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Abstract: We have studied the nature of classical work (Wc) and thermodynamic 
work (W) fluctuations in systems driven out of equilibrium both in transient and time 
periodic steady state. As the observation time of trajectory increases, we show that 
the number of trajectories which exhibit excursions away from the typical behaviour 
i.e., Wc < 0, < AF and dissipated heat Q < decreases as anticipated for 
macroscopic time scales. Analytical expressions for such trajectories are obtained. 
Trajectory for which Wc < may not correspond to W < AF or Q < 0. The 
applicability of steady state fluctuation theorems are discussed in our linear as well 
as nonlinear models. 
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I. INTRODUCTION 



Over last decade, nonequilibrium fluctuation theorems have attracted much interest. 
They provide relations for physical quantities such as work, heat and entropy in driven 
nonequilibrium systems, independent of the nature of driving [1-17]. The fluctuation 
relations are statements about the symmetry of the distributions of the physical quan- 
tities around zero and not around the maximum. These distribution functions exhibit 
flnite weight for negative values for physical quantities of interest, which are usually rare 
and are related to transient second law violating contributions at microscopic length 
scales. Fluctuation theorems quantify the probability of those non-equilibrium trajec- 
tories, taken individually, violate some of the inequalities of the thermodynamics. For 
average quantities, these theorems lead to inequalities consistent with the second law of 
thermodynamics. In our present study, we are mainly concerned with Jarzynski Equal- 
ity (JE) [10,11] for thermodynamic work (W) and Bochkov-Kuzovlev identity (BKI) for 
classical work (Wc) [14-17]. JE relates the nonequilibrium work done with equilibrium 
free energy difference (AF — Fb — Fa) between two thermodynamic states. These two 
thermodynamic states are uniquely characterized by the initial (A) and flnal (B) values 
of the time dependent protocol. JE is given by 

(e-/^^) = e-^^^, (1) 

where W is the thermodynamic work. The angular bracket (• • •) denotes the average 
over an ensemble of realizations, starting from initial equilibrium conflgurations. W 
is deflned as W — Jq ^^^'^^ dt, where U{x, t) is the effective potential of the system. 
Incidentally, the JE flnds its counterpart in the relation known as Bochkov-Kuzovlev 
identity which has been proposed much earlier [14-17] 

(e-^^=) = 1 (2) 

where Wc is the classical work. If the potential is decomposed as U{x,t) = 
Uo{x) + Up{x,t), the classical work done over a time interval r is deflned as Wc = 



— Jq — ^^xdt. For the validity of BKI, the initial equilibrium distribution must 

'^o(^o) 

correspond to the time independent potential {Pe{xo) = Ne ''b^ , being the nor- 
malization constant). Seifert has shown that both JE and BKI are special cases of 
more general result, obtained within a framework that defines entropy production in 
nonequihbrium state [18,19]. When the internal energy of the system is defined by the 
bare Hamiltonian/potential, then classical work is defined as the work performed by the 
application of an external force that affects the system in a fixed landscape of the bare 
Hamiltonian. To define thermodynamic work, external time dependent perturbation is 
treated as the time dependent contribution to the internal energy of the system [20]. 
The work done by the system on the external bodies which produce a change in the 
Hamiltonian is known as thermodynamic work W. It turns out that W is more useful 
than Wc for the reconstruction of free energy landscapes. Using Jensen's inequality, 
from eqns. (1) and (2) it follows that (W) > AF, (Wc) > 0, which are consistent with 
the physics at macroscopic scales. Observations of realizations wherein W < AF or 
Wc < arc treated as excursions away from the typical macroscopic behaviour. The 
finite time trajectories for which W < AF, are sometimes referred to as transient second 
law violating trajectories [21]. In our work, we have also studied decrease in the number 
of such trajectories as a function of time of observations and some pertinent questions 
are raised. Nature of the fluctuations in W and Wc are studied in a nonlinear system 
exhibiting stochastic resonance (SR) [22-29]. Finally we have analyzed the vahdity of 
SSFT for W and Wc- 

II. DRIVEN LINEAR SYSTEMS 
A. Model 

We consider a linear model, i.e., a overdamped Brownian particle in a one dimensional 
harmonic potential Uo{x) — \kx'^ in the presence of an external time dependent 



potential Up — —Ax{t) Sin.{ujt + 0), where A , u) are the amphtude, frequency of the 
external drive respectively and (f) is an initial phase. The dynamics of the particle is 
described by the Langevin's equation 



where U{x,t) — Uq{x) + Up{x,t). The random force field ^(t) is a zero mean Gaussian 
white noise, i.e., {^{t)^{t )) = 2D6{t -t),D = -fkeT is the noise strength of the 
medium. Here 7 is the friction coefficient, ks is the Boltzmann constant and T is the 
absolute temperature of the bath. In the following we use a dimensionless form of 
eqn.(3), namely 

|-^-«')- w 

All the calculated quantities are in dimensionless form. 



B. Thermodynamic work distributions: 

Thermodynamic work(Vr) done on the system by the external drive over a time t is 
given 

W = -Auj / xit) cosiiot + (t))dt . (5) 
Jo 

The formal solution of eqn.(4) is 

xit) = XQe-icpi-kt) + / dt ex.p(-k(t - t))[Asm(u;t + 0) + ^(t)]. (6) 

Jo 

Where Xq is the initial coordinate of the particle. The initial distribution for xq is 
assumed to be the equilibrium canonical distribution, 

P^(xo) = ^/^ex.p{-!f{xo -jsm{(j))f). Prom eqns.(5) and (6), it follows that 
thermodynamic work done is a linear functional of the Gaussian variable ^ (t) . Hence 



the distribution P(W) of work W is a, Gaussian [31-33]. By a simple algebra one can 
evaluate P{W) analytically. The full probabihty distribution P{W) is given by 



{W-{W)f 
2^^ 



(7) 



where the analytical expressions for the mean thermodynamic work done{{W)) is given 
in Appendix-A. The variance cr^ = {W'^) — {W)"^ is given by 

(W) - {Wr = ^{W) + ^sm\u;t + <f>). (8) 

The thermodynamic work done over a time interval r satisfies JE, 

{e-^W) = e-^^^ (9) 

where AF = — ^ sin^(ci;r + 0), is the free energy difference between two 
thermodynamic states with potentials U{xo) = ^kx^ — Axosm(4') and 
U {x, t) = ^kx^ — Ax{t) sin(a;T + 0) respectively . 

We furthur study the statistics of the thermodynamic work done, W, in the time 
asymptotic regime (i.e, in the limit t — > oo). In this limit, probability distributions are 
time periodic with a period of magnitude The average work done((Wr)) over a 
time of observation r in the time periodic steady state is 

(Wr) = }im[{W{t + T))-{Wm 

t— >oo 

(10) 



2(F + cx;2)' 

The variance((T^) of W-,- averaged over a time of observation r is given by 

^^{W^) + A{t) (11) 

The analytical expression for A(r) is given in Appendix-C. The distribution for Wr is 
Gaussian. However, it does not satisfy the fluctuation dissipation relation, namely 



al — [4]. Hence SSFT is not valid over small time of observation. However, 

(Wr) is an extensive quantity in r, whereas A(t) saturates to a finite value for t » |. 
For large r, ^{Wr) ^ ^(t) and the contribution of A(r) to the variance can be 
ignored in eqn.(ll). Under this approximation SSFT holds. Here SSFT implies 



= e^""^ (12) 



P{-W,) 

C. Classical work distributions: 

Classical work {Wc) done on the system by the external drive As\n{ijjt + 0) over a 
duration of time t is given by [14,15] 

Wc ^ I Ashiiijot + 4>)x{t)dt 
Jo 

^ W + [A sm{ujt + ^)x{t) - A sin((/))a;(0)] (13) 

Note that classical work (Wc) differs from the thermodynamic work(W) by a boundary 
term. From eqns.(13) and (6) it follows that Wc is a hnear functional of a Gaussian 
variable ^(t), consequently P{Wc) is a Gaussian and is given by 



PiWc) = , exp 



iw.-{Wc)r ^^^^ 



c 



where (VFc) is the mean classical work done over a time interval t and cr^ is the 
fluctuation or variance of Wc- These quantities can be evaluated analytically. The 
expression for average or mean classical work((VFc)) is given in Appendix-B. 
The variance in classical work (cr^) over a transient time t is given by 



al = {W'c) - {Wcf 



l^Wc) (15) 



To calculate Wc and a^, we assumed initially the system to be in equilibrium with 
distribution P^{xo) = \f{^)e'^^^- The expression for P{Wc) is consistent with BKI. 

(e-/^^=) = 1. (16) 

The average classical work ((Wc-)) done over a time of observation r in the time 
periodic steady state, is given by 

which has the same magnitude as that of {Wr) and is independent of (f). The 
distribution P{Wct) is Gaussian. The variance {(Jct)) is given by 

(Tcr=|(We.)+Ai(T) (18) 

An expression for Ai(r) is given in Appendix-D. For r ^ |, Ai(r) saturates to a finite 
value. For large r, for which ^{Wct) ^ Ai(r) one can ignore Ai(r) in eqn.(18). Under 
this condition the classical work satisfies SSFT. 



In fig(l) and (2), we have plotted the probability distributions of W and Wc for various 
values of time periods in the transient regime using our analytical results 
(eqns.(7),(14)). Other parameters are given in figure captions. For short time 
t — 0.25Tp {Tp being the time period of magnitude ^), we notice that most of the 
weight of W is located in the negative side. The weight towards negative side for 
values W < AF comes from the so called transient second law violating trajectories 
[21]. As we increase the time of observations, the total weight for W < AF (i.e., area 
under the curve P{W) between W — — oo and W — AF) decreases (see fig(l)). We 
would like to emphasize that depending on the protocol and potential, one may obtain 
peak in P{W) at values W < AF, i.e., most probable value of W is inconsistent with 
the classical macroscopic or thermodynamic results. However, (W) > AF is always 
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FIG. 1: The distribution P{W) as a function of transient time t (in units of Tp) for parameters 
A = 0.3, D = 0.1, u) = 0.1 and T„ = ^. 

satisfied. In contrast, Wc does not exhibit a peak in the negative region. We have 
checked this separately for different systems for different protocols. 
The physical quantities such as work (W), heat (Q), total entropy (AStot) and internal 
energy (At/) can be calculated using the method of stochastic energetics [34]. For 
details, we refer to [26,27,28,35]. In fig(3) we have plotted number of trajectories which 
do not satisfy our conventional wisdom at macroscopic scale, namely W < AF, Wc < 
and heat dissipated, Q < 0. For each individual trajectory, W, Wc and Q are 
calculated. In our simulation, we have generated 10^ trajectories. Number of 
trajectories for which Wc < 0, W < AF and Q < are all different. Trajectory for 
which Wc < need not correspond to < AF or Q < 0. It can be readily shown that 
if is the total number of observed trajectories for Wc, then the number of trajectories 



w 



FIG. 2: The distribution P{Wc) as a function of transient time t (in units of Tp) for parameters 



A = 0.3, D = 0.1, cj = 0.1 and Tp = ^. 



A^' for which Wc < are given by A^' = YErfc( [ ). In the time asymptotic regime 



{Wc) scales with t and the decay of A^' with time t is given by A^' = ^ exp(— ct), where 
a and c are constants. For our present problem a = J ^'^4^2 " ^ and — — Ajl^ 



The number of trajectories for which W < AF in the large time limit, decay with the 
same functional form as that for Wc < 0. Hence there is no correlation amongst 
trajectories in regard to the transient violations in W, Wc and Q. For our simple linear 
problem and protocol (up to time periods ~ 12), number of trajectories for which 
Q < are greater than W < AF which in turn are greater than Wc < 0. This is not a 
general rule. Depending on the system and protocol, different possibilities exist. 
However, one common observation is that the number of trajectories which defy our 
general classical notion for W, Wc and Q decrease monotonically as time of observation 
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FIG. 3: Number of trajectories for W < 0, Q < and Wd < as a function of time t (in 
units of Tp). The physical parameters are same as in fig(l) and (2). 



is increased, leading to classical thermodynamic behaviour in the time asymptotic 
regime as is evident from fig(3). To put the probability distributions of physical 
quantities W, Wc, Q and AStot in the same perspective we have generated 10^ 
trajectories of Brownian particle. Total entropy production AStot comprises of two 
parts namely system entropy and medium entropy. For details we refer to [18,19]. For 
each stochastic trajectory we calculate W, Wc, Q and AStot and using these values, the 
obtained P{W), P{Wc), P{Q) and P{AStot) are plotted in fig (4). Areas under the 
curves for W < AF, Wc < 0, Q < and AStot < are different. In our linear problem 
the distribution for P{W), P{Wc) and P{AStot) are Gaussian, while that for P{Q) is 
non Gaussian. It may be noted that by appropriately choosing the system and 
protocol most probable value of P{AStot) can be shifted to the negative side, yet 
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FIG. 4: The distributions P{AStot), P{W), P{Wc) and P{Q) for the parameters A = 0.3, 



D = 0.1, uj = 0.1 and for t = ^. 



{AStot) > 0. However, peak in P{Q) corresponding to the most probable value occurs 
at positive values of Q only. Similar behaviour in regard to P{W) and P{Wc) are 
pointed earlier. In the next section, we study the nature of fluctuations in W and Wc 
for a nonlinear driven system in the time periodic steady state. 



III. DRIVEN NONLINEAR SYSTEMS 



To this end we consider the motion of a Brownian particle in a bistable system under 
the action of an external ac force. The total potential is given by 
U{x,t) = Uo{x) + Up{x,t). The static double well potential is Uo{ and 
the time dependent potential is Up{x,t) = —Ax(t) cos(u;t + (/)). This system is shown 
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FIG. 5: The distribution P{Ws) over a single period for various values of initial phase 4>. The 
other parameters are A = 0.1, D = 0.12 and uj = 0.1. Inset shows plot of (Wg) vs D for 
different value of (p for the same physical parameters. 



to exhibit the well known phenomenon of stochastic resonance (SR) [22-28]. 
Theoretical study on thermodynamic work and heat have been shown to satisfy SSFT 
[26,27]. Experiments in connection with fluctuation theorems have been carried out in 
systems exhibiting SR [28,29]. We have numerically evaluated the classical work and 
thermodynamic work distributions over a single period/large number of periods. For 
description of numerical method we refer to [26,27]. 

In fig(5) and (6) we have plotted probability distributions of thermodynamic work Wg 
and classical work Wcs evaluated over a single period in the time periodic regime for 
various values of initial phases. For a given fixed phase, distributions P{Ws) and 
P{Wcs) are different. Only for phase <p = ^, both distributions are identical. For this 




FIG. 6: The distribution P{Wcs) over a single period for various values of initial phase (p. The 
other parameters are A = 0.1, D = 0.12 and uj = 0.1. Inset shows plot of {Wcs) vs D for 
different values of (p for the same physical parameters. 



particular phase, it can be readily seen that thermodynamic work and classical work 
are identical when evaluated over a period(s). The multipeaked structure in the 
distributions arise from the interwell and intrawell dynamics of the particle in a double 
well system. For details see Ref. [27,28]. Unlike the sensitivity of P(Wcs) and P(Ws), 
{Wg) and (Wcs) over a cycle do not depend on the initial phase and are identical, i.e., 
(Wcs) = i^s)- They are plotted in the inset of fig(5) and (6). (Ws)(= {Wcs)) exhibits a 
peak as a function of noise strength D. This phenomenon is referred to as SR 
[25,26,27]. At the value of noise strength corresponding to the peak value in (Wc), the 
random hops of the Brownian particle between the two wells get synchronized with the 
external drive[22,28]. 
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FIG. 7: The distributions P{Wn) for different periods (n) in tlie time periodic steady state. 
In the inset P{Wn) is plotted along with Gaussian fit f{W) for n = 15. From the Gaussian 

2 

fit f(W), the fluctuation ratio has been calculated as ■> f^,, , = 1.022. Here D = 0.12, A = 0.1 
and uj = 0.1. 



In fig(7) and (8), we have plotted P{Wn) and P(Wcn) versus thermodynamic work 
done over n cycles, Wn and classical work done over n cycles. Wen respectively. In both 
these cases, with increasing n, multipeaked distributions become smoother and for 
large n, they both evolve towards Gaussian distribution. Corresponding Gaussian fit 
for P{Wn) and P{Wcn) for n = 15 are shown in the inset of fig(7) and (8) respectively. 
For P{Wn) and for n = 15, from the Gaussian fit we get variance o"^ = 0.3206, 
(Wn) = 1.30701 and the corresponding fluctuation ratio, -r^rr = 1.022, which is close 
to unity (within our numerical accuracy). The Gaussian nature of Wen along with 
fluctuation ratio being 1 implies validity of SSFT. Similar conclusions can be made for 
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FIG. 8: The distributions P{Wcn) for n periods versus Wen in the time periodic steady state. 
Inset shows the plot of P{Wcn) along with it's Gaussian fit f{Wc) for n = 15. From the 



Gaussian fit f{Wc), the fluctuation ratio has been calculated as 
parameters are same as in fig(7). 



1.005. The other 



P{Wcn) for n = 15, where the Gaussian fit gives cr^ = 0.31525, and hence the 
fluctuation ratio (1.005) is close to unity. The presence of non-Gaussian tails at large 
value of Wen and Wn are not ruled out. However, numerically it is difficult to detect 
them. It may also be noted that time interval after which SSFT is obeyed for W and 
Wc may be different and depends of the physical parameters. 



IV. CONCLUSION 



We have studied fluctuations in W and Wc analytically for both transient and time 
periodic steady states in case of a linear model. We have also shown that as the 
observation time of trajectories increases, the number of trajectories which exhibit 
atypical behaviour (namely Wc and Q being negative, W being less than AF) 
decreases. Thus in the limit of large time of observations, macroscopic thermodynamic 
behaviour results for physical quantities. We have discussed the validity of SSFT for 
both W and Wc in our linear as well as nonhnear models. 
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VI. APPENDIX-A 



(W) = 



4(A;2+u;2) 



[cos{2{iot + 0)) - cos(20)] 



+ /, n — — exp{—kt)[—k cos{ut + 0) + u!sm{u)t + 0)](A;sin(0) — u cos(0)) 



+ 



A^u;sin((/)) 



(A;sin(0) — a;cos(0))(A;cos(0) — a;sin(0)) 
[exp(— A;t){— A;cos(6i;t + 0) 



A;2 + cu2 

+ uj sm{ut + 4>)} + {k cos{(f)) — u sm{(f)))] 



(19) 



VII. APPENDIX-B 



^^^^ ^ i(^^^W2M + <^))-cos(2<^)] 

+ 77q ov7 exp(— A;t)[— A; cos(c<;t + 0) + a;sin(c<;t + 4>)](ksm(4>) — u cos(0)) 

+ (fc2 ^^2)2 (^^^^('^) -^cos((/)))(/ccos(0) -wsin(0)) 

+ -ex.p(—kt)sm(u;t + (j))[u;cos((l)) — ksm((f))]. (20) 



VIII. APPENDIX-C 

The expression for A(r) is given by 

942 2 

A(r) = -(l^.) - ^^^^^^-_^(A:2cos2(0) -a;^sin2(</)))(l -exp(-/cr)) (21) 



IX. APPENDIX-D 



The analytical expression for Ai(r) is given by 
2 2A'uj ^ 



Ai(r) = -{W,r) + Tl^rnr—^ik\os\(P)-u^sm\4>)){l + exp{-kT)) 



AA^uj^k 2A^ . 2, , 2^2 (P- a; 



2N 



+ -— sin^(0) + ^wTTo ^ sin^(0) exp(-A;-r) (22) 



l3{K^ + uj^Y k/3 '"^^ kp{k^ + u;^ 
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